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ABSTRACT 

From July 2008 to October 2009, the Gamma-ray Bm-st Monitor (GBM) on board the Fermi 
Gamma-ray Space Telescope (FGST) has detected 320 Gamma- Ray Bursts (GRBs). About 20% 
of these events are classified as short based on their Tgo duration below 2 s. We present here for 
the first time time-resolved spectroscopy at timescales as short as 2 ms for the three brightest short 
GRBs observed with GBM. The time-integrated spectra of the events deviate from the Band function, 
indicating the existence of an additional spectral component, which can be fit by a power-law with 
index ^ —1.5. The time-integrated E'peak values exceed 2 MeV for two of the bursts, and are well 
above the values observed in the brightest long GRBs. Their -Epeak values and their low-energy power- 
law indices (a) confirm that short GRBs are harder than long ones. We find that short GRBs are 
very similar to long ones, but with light curves contracted in time and with harder spectra stretched 
towards higher energies. In our time-resolved spectroscopy analysis, we find that the -Epeak values 
range from a few tens of keV up to more than 6 MeV. In general, the hardness evolutions during the 
bursts follows their fiux/intensity variations, similar to long bursts. However, we do not always see 
the Epcak leading the light-curve rises, and we confirm the zero/short average light-curve spectral lag 
below 1 MeV, already established for short GRBs. We also find that the time-resolved low-energy 
power-law indices of the Band function mostly violate the limits imposed by the synchrotron models 
for both slow and fast electron cooling and may require additional emission processes to explain the 
data. Finally, we interpreted these observations in the context of the current existing models and 
emission mechanisms for the prompt emission of GRBs. 

Subject headings: gamma-ray burst: general, gamma rays: general, radiation mechanisms: non- 
thermal 



1. INTRODUCTION 

The bimodal distribution of Gamma-Ray Burst (GRB) 
durations, with a mi nimum around 2 s, wa s first 
clearly established by iKouveliotou et all (|1993 ^ using 
data from the Burst And Transient Source Experiment 
(BATSE) on-board the Compton Gamma- Ray Observa- 
tory [CGRO). Further evidence for two separate popu- 
lations was a tendency identified in the same data for 
the short bursts to have harder spectra than the long 
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bursts (IKouveliotou et all 119931: iPaciesas et all 120031) . 
A natural consequence of these observations was the 
speculation that the separate populations might have 
distinct physical origins. Since then, detailed multi- 
wavelength observations of GRB afterglows have sup- 
ported a physical picture in which long GRBs are 
produced by the gravitatio nal collapse of very mas- 
sive stars C Wooslcy 1993: Ma cFadven and Wooslevlll999l: 
IWooslev and Heger. .2006. ) and short GRBs are the re- 
sult of mergers of a neutron st ar with anoth e r neutron 
star or a bla ck hole IPaczvnskH (BH - I1986D ; IRosswod 
(BH - I2OOI : alternate sce narios include the merger of 
a BH wit h a white dwarf (iFrver et al.l |1999() or a he- 
lium star ()Frver and Wooslevll998() . Finallv. short GRBs 
could a lso be produced by a variant of the collapsar 
model (| Zhang et al.l (I2003f) - for ext ens ive reviews o n 
long and short GRBs see lNakari (|2007D or lPiranI (|2004D ). 

Additional observational evidence differentiating the 
burst classes came from spectral lags, estimated from 
cross-correlations of burst time-profiles between differ- 
ent energy bands within the 25-1000 keV BATSEl\jSi^ e 
Area Detectors (LADs) energy range (|Norris et al.ll200Q[ ) . 
The distributions of spectral lags of short and long GRBs 
diff er significantly. For long GRBs the typi cal lag is ^100 
ms (jNorris et al.ll2002t iHakkila et al.ll2007t ) with the hard 
photons arriving earlier (hard-to-soft evolution), whereas 
short GRBs typically have lag values consistent with 
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zero (|Norris and Bonnel]|[200llYi et al.ll2006D . Interest- 
ingly, however, recent results with the Fermi Gamma Ray 
Space Telescope {FGST) indicate that in the GeV range, 
photons always arrive later than the prompt gamma-rays 
(10 — 1000 keV), thus reversing the trend to a soft-to-hard 
evolut io n both in long and in short GRBs (jAbdo and alJ 
I2009cllbl : lAckermann and al.ll2010al) . 

Whether the nature of the central engine differs be- 
tween the two classes is not clear. However, the general 
picture of a collimated ultra-relativistic outflow seems 
to apply in both cases and theories for generating the 
prompt and afterglow e mission gen erally do not dis- 
tinguish between them (|Piranl [2004D . Although time- 
resolved spectral analysis may provide significant the- 
oretical constraints, most short GRBs are too weak 
for detailed spectral analysis and one is usually re- 
stricted to simplified parametrizations such as hardness 
ratios and spectral lags, which provide rather crude con- 
straints for theoretical models. Prompt GRB spectra can 
gener ally be fit a dequately with the Band GRB func- 
tion ([Band et al.l |T993i) , which is basically two power 
laws joined smoothly at a break energy, -Epeak, indicat- 
ing the maximum of the vF^, spectrum. However, the 
high-energy power-law is often not well measured or con- 
strained, especially for weak and/or short bursts, and a 
single power-law with a high - energy cutoff can be used in- 
stead. [Pacieiai^^F^I] (|2003f) investigated the distribution 
of the spectral parameters of BATSE GRBs versus dura- 
tion and found that short GRBs differ from long GRBs by 
having both high er E'pcak values a.nd harder low -energy 
power- law indices. iGhirlanda et al.l (|2004l2009bl ) studied 
a sample of relatively bright BATSE GRBs and found 
that the hardness difference was attributable mainly 
to the low-energy power-law index with no difference 
in -Epcak values. They also compared the short-GRB 
spectral-parameters with the first few seconds of long 
GRBs and found no significant difference. 

The Gamma-ray Burst Monitor (GBM) on board 
FGST has been detecting GRBs since July 2008. From 
July 2008 to October 2009, GBM triggered on about 
320 GRBs, of which approximately 20% have durations 
shorter than 2 s (Paciesas et al., in preparation). With 
the unprecedented time resolution of GBM, down to 
about 2 /is, we have the ability for the first time to 
perform fine-time resolved spectroscopy of the bright- 
est short GRBs at a millisecond timescale. Such de- 
tailed analysis is a key point in understanding the differ- 
ences and similarities in the behavior of the two classes 
of GRBs and their emission mechanisms. 

In this paper we have selected a subset of the three 
brightest short events detected with GBM and performed 
detailed temporal and spectral analyses at millisecond 
timescales. In section [5J we describe the instrument, 
pointing out its assets in the context of the following 
analysis, as well as the data set we used. Section [3] is de- 
voted to the timing analysis with time-lag and duration 
studies. The time-integrated and time-resolved spectral 
analysis of the three selected events are presented in sec- 
tion Section [S] is dedicated to the theoretical interpre- 
tation of our observations. 

2. INSTRUMENTATION AND DATA 

FGST w'ds launched on 2008 June 11 with two instru- 
ments on board: the Large Area Telescope (LAT) and 



GBM. The latter instrument is an array of scintillation 
detectors sensitive to gamma-rays in the energy range 8 
keV to 40 MeV. Twelve sodium iodide (Nal) detectors, 
12.7 cm in diameter by 1.27 cm thick, cover energies up 
to 1 MeV and are used to determine burst locations. Two 
bismuth germanate (BGO) detectors, 12.7 cm diameter 
by 12.7 cm thick, cover energies above 200 keV. The high 
sensitivity of GBM above 1 MeV is ideal for the study 
of hard GRBs. Moreover, the broad energy range cov- 
ered by the combination of the Nal and BGO detectors 
allows a good measurement of the GRB prompt spectral 
parameters over a wide range. The broad energy range 
of GBM also permits investigation of possible deviations 
of the GRB spectra from standard empirical models (e.g. 
the Band function), which were in the past applied within 
a narrower energy domain . The instrument is described 
in detail bv lMeegan et al.l ([2009). 

The very high temporal resolution of GBM is a major 
asset for the study of short events. The GBM onboard 
software incorporates burst triggering on time scales as 
short as 16 ms. All triggers generate Time- Tagged Event 
data (TTE) consisting of arrival times and the deposited 
energy for individual photons/events from each of the 14 
detectors with a temporal resolution of 2 fis. A pre- 
burst ring buffer saves 512000 events before the trig- 
ger, which corresponds to ~30 seconds at background 
rates. TTE data are produced for 300 seconds after the 
trigger, thus all short bursts have full temporal cover- 
age. The energy is digitized into 128 channels, pseudo- 
logarithmically spaced to provide channel widths less 
than the detector resolution up to 12 MeV, though TTE 
data are available at cruder spectral resolution up to 40 
MeV. 

Each of the GBM detectors has a fixed dead time per 
event of 2.6 /iS, which is corrected in the present study. 
In addition, the spacecraft telemetry bus imposes a maxi- 
mum total event rate (sum of all detectors) of 375k events 
per second. Neither of these limitations is significant for 
the bursts described here. 

GBM detected about 70 short bursts between 2008, 
July 14 and 2009, October 1. Although its broad en- 
ergy coverage, its good sensitivity above 1 MeV (due to 
the thickness of the BGO detectors), and its fine spec- 
tral resolution make GBM an ideal instrument to study 
GRB spectra, the instrument has a lower effective area 
than the BATSE/L ADs and a higher background than 
the 5'wii/t/Burst Alert Telescope (BAT); as a result most 
short bursts are not fluent enough for spectral analysis to 
be very constraining, particularly if one wishes to study 
spectral evolution over time. This low fluence problem is 
aggravated by their typically higher -Bpcak values, which 
make both -Epoak and the high-energy index of the Band 
function (/3) more difficult to measure for short bursts. In 
the following study, we have selected the GRBs detected 
with GBM from the start of the mission untill October 
2009 with a fiuence above 2 x 10~^ ergs cm~^ between 8 
and 1000 keV. This fluence threshold assures the possi- 
bility of performing a time-resolved spectral analysis at 
the milli-second time-scale during the entire main emis- 
sion of the burst. This selection criterion results in a 
sample of 3 bursts: GRBs 090227B, 090228 and 090510. 
Time-resolved analyses are still possible with other short 
GRBs detected with GBM, albeit in a coarser temporal 
and spectral resolution. The fourth brightest short GRB 
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during this time period, GRB 081216, with 1.4 x IQ-^ 
erg cm~^, has 2.8 times less fluence than GRB 090510, 
which limits the temporal resolution at which spectral 
analysis can be performed during the weaker portions of 
the GRB, requiring wider intervals that merge peaks and 
valleys of the light curve. 

GRB 090510 triggered G EM and LAT on 2009 
May 10 at 00:22 :59.97 UT (iGuiriec et al. 
00:23:01.22 UT (p^^l [200l lOmodei et a" 
spectively. The burst was also detected with 
Swift/BAT (jHoversten et all 120091 ) and located with 
the Swift/mtTa Violet O ptical Telescope (UVOT - 
iKuin and Hoverstenl I2OO90 precisely enough to perform 
follow-up observations with optical ground-based tele- 
scopes. The Very Large Telescope/FOcal Reducer and 
low dispersion Spectrograph (VLT/F0 RS2) measured 
a spectroscopic rcdshift of 0.903±0.003 ([McBreen et al.l 
\20m IRau et al.. 2009) . 

GB M detected consecuti vely GRB q90227B (|Guiried 
[20091) and GRB 090228 (Ivon KienhnI [20091) on 2009 
February 27 at 18:31:01.41 UT and on 2009 February 28 
at 04:53:20.91 UT, respectively. These bright events were 
also seen by several other instruments: Konus-Wind, 
Suzaku-WAM, INTEGRAL SPLACS,and MESSEN- 
GER for GRB 090227B; Konus-Wind, Mars Odyssey, 
MESSENGER and Agile-MCAL for GRB 090228. This 
allowed reconstruction of the arrival direction of their 
photons by the Inter-Planetary Network (IPN) and pro- 
vided an error box for the position of each GRB (Hurley 
et al.: 3'''^ IPN Catalog in preparation). These loca- 
tions typically achieve more precision than is possible 
with GBM alone. The best locations for all three GRBs 
are reported in Table [1] The set of the brightest GBM 
detectors with an angle to each source below 80° and a 
source view not blocked by the spacecraft's components 
(solar panels, radiators, LAT, etc.) are reported in Ta- 
bled 



TABLE 1 

Best locations obtained for GRBs 090227B, 090228 and 
090510. 



GRB 


Instrumcnt(s) 


RA (°) 


DEC (°) 


ERR 


GRB 090227B 
GRB 090228 
GRB 090510 


IPN 
IPN 

Swift/VYOT 


15.6 
98.3 
333.6 


26.4 
-28.4 
-26.6 


0.6" 
O-OOl?"^ 
1.5' 




" in square degrees; in arc seconds 



TABLE 2 

Detectors used for the spectral analysis in each GRB. 
We used the brightest detectors with an angle to the 

SOURCE below 80° FOR THE NaI DETECTORS AND THE CLOSEST 
DETECTOR TO THE SOURCE FOR THE BGOs. 



GRB 


Nal detectors 


BGO detectors 


GRB 090227B 
GRB 090228 
GRB 090510 


nO, nl, n2, n5 
nO, nl, n2, nS, n5 
n3, n6, n7, n8, n9 


bO 
bO 
bO, bl 



3. OBSERVATIONS: I. TEMPORAL ANALYSIS 
3.1. Light-Curves and Time lags 

The high temporal resolution (2 ms) light-curves of all 
three GRBs are shown in figures [T] to [3] in two energy 
ranges, 8 — 200 keV {panels a) and 1-38 MeV {panels h) , 
respectively. GRB 090510 (Figure (H]) triggered GBM at 
To on a very soft and weak pulse, not shown in the plot, 
which was followed by an interval of about 0.5 s without 
significant emission. Figure [3| shows the main emission 
from this burst starting at Tq + 0.5 s. 

For all three GRBs, the higher energy light curves 
exhibit much sharper structures than the lower energy 
curves, where the emission is smoother. GRBs 090227B 
and 090510 show very fast variability with multiple spikes 
in all energy bands; GRB 090228 is much simpler with 
only two main peaks. Overall, the signal rises much faster 
than it decays. A visual comparison of the high and low 
energy light curves reveals a trend of soft-hard-soft evolu- 
tion in all events, which is also supported by our spectral 
fits as discussed in section B. 2. II The most dramatic evo- 
lution is demonstrated in the data of GRB 090510 (Fig- 
ure [3l). 
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Fig. 1. — Light curves of GRB 090227B in two energy bands 
{Panel a: 8 keV to 200 keV, Nal detectors) and (Panel b: 1 
MeV to 38 MeV, BGO detectors) with 2 ms time resolution. The 
count rates are background substracted. Two bottom panels: The 
same light-curves with variable time bins (histograms), optimized 
for time-resolved spectroscopy. The Band function peak energy, 
^'pcaki is plotted over the lightcurve for each time interval. 



We have estimated the time averaged spectral lags for 
each of the events using a cro ss-correlation technique 
(jCheng et al.l Il995t iBandl I1997D . We use the Pearson 
cross-correlation function to estimate the degree of cor- 
relation of two independent time series. For each detector 
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Fig. 2. — Light curves of GRB 090228 in two energy bands {Panel 
a: 8 keV to 200 keV, Nal detectors) and (Panel b: 1 MeV to 38 
MeV, BGO detectors) with 2 ms time resolution. The count rates 
are background substracted. Two bottom panels: The same light- 
curves with variable time bins (histograms), optimized for time- 
resolved spectroscopy. The Band function peak energy, i?pcaki is 
plotted over the lightcurve for each time interval. 
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Fig. 3. — Light curves of GRB 090510 in two energy bands (Panel 
a: 8 keV to 200 keV, Nal detectors) and (Panel b: 1 MeV to 38 
MeV, BGO detectors) with 2 ms time resolution. The count rates 
are background substracted. Two bottom panels: The same light- 
curves with variable time bins (histograms), optimized for time- 
resolved spectroscopy. The Band function peak energy, Speak > is 
plotted over the lightcurve for each time interval. 



type (Nal and BGO), we generated the GRB light curves 
using the TTE data divided into 6 energy bands: 4-20 
keV, 20-41 keV, 41-70 keV, 70-212 keV, 212-510 keV, 
and 510-985 keV (Nals); 109-212 keV, 212-517 keV, 
517-996 keV, 996-2048 keV, 2.0-4.0 MeV, and 4-45.5 
MeV (BGOs). We then measured the lags between the 
first Nal energy band and each of the others. To evaluate 
the errors on the lags, we performed the same analysis 
with simulated light-curves generated from the real light 
curves by adding Poisson fluctuations based on the count 
rate in each time bin. The lag errors are estimated from 
the distribution of 100 such trials. The results are pre- 
sented in Figure m Negative (positive) numbers indicate 
that the higher energy photons lead (lag) the low energy 
photons. 

In general, short GRBs are expec ted to exhibit small 
or no spectral lags below 1-2 MeV (|Norris and Bonnelll 
I2006f l. Figure H shows that for GRB 090227B there is a 
constant lag of ~ 15 ms {l-2a level) between the higher 
energy (30—3000 keV) photons and those at ~8-20 keV 
. For GRB 090228, although the lag values are all pos- 
itive (indicating a lag of higher energies with respect to 
low energies) the errors are too large to allow a signifi- 
cant measurement. Finally, when the burst has sufficient 
counts in the higher energy bands (> few MeV), for 
GRB 090510, we measure a positive lag between 8 keV 
to 3 MeV and the highest energy band (> 3 MeV). This 
result is also supported by the higher energy LAT data, 
as already pub lished in the supplementary information 
of lAbdo andal.. (,2009a) . The energy range 8-20 keV is 



outside that used in iNorris et al.l ()2000t ). The canonical 
energy channels used to provide a discriminant between 
short and long bursts show no eviden ce for lags in any of 
the th ree bursts, consistent with the INorris and Bonnelll 
(|2006l ) observation that short bursts show small/zero- 
lags. 

The average lags measured above indicate an overall 
good alignment of the light-curves from 20 keV to sev- 
eral MeV. However, local time lags and leads could still 
be present between individual pulses, as pr eviously ob- 
served for a iiumbe r of long BATSE GRBs (|RvdeH2005l : 
iHakkila et al.ll2008l) . In particular, we note that, as also 
seen in the two upper panels (a and b) of Figure [21 in 
GRB 090228 the peak of the brightest pulse in 8-200 keV 
lags the peak in the >1 MeV band. Such local delays 
cancel out when integrated over the entire light curve, 
resulting in a small or zero lag. 

3.2. Durations 

We have performed a detailed temporal analysis of 
the three selected GRB s, app lying the technique de- 
scribed in iKoshut et al.l (|1996[ ) to determine the Tgo 
and T50 durations of all three events in the 50-300 
keV energy range using the TTE data type (see sec- 
tion [2|). Tgo (T50) is the time between accumulating 
5% and 95% (2 5% and 75%) of the co unts associated 
with the GRB (jKouveliotou et al.lll993f) . We describe 
below only the results of this analysis for GRBs 090227B 
and 090228; the duration analysis of GRB 090510 has 
been described in the Supplementary Material of the pa- 
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Fig. 4.— Spectral lags for GRB 090227B (top panel), 
GRB 090228 (middle panel) and GRB 090510 (bottom panel). 
All lags were computed with respect to the 8-20 keV Nal energy 
band (see text for details). The temporal resolution was 2 ms for 
GRBs 090227B and 090228, and 1 ms for GRB 090510. 



per bv lAbdo and al.l (|2009a^. For both events we used in 
our analysis the combined light curves of the three Nal 
detectors with the smallest zenith angles to the source. 
For GRB 090227B these are Nal 0, 1, and 2 with angles 
ranging from 19° to 52°; for GRB 090228 these are Nal 
0, 1, and 3, ranging from 8° to 42°. Table [3] contains 
Tgo and T50 results for all three events. 





Burst 


Tbo (s) 


Tso (s) 




GRB 090227B 


2.38 


0.09 




GRB 090228 


0.34 


0.05 




GRB 090510 


2.1 


0.2 



TABLE 3 

Tgo AND T50 DURATIONS. 



Further, we computed the fluence hardness ratios (HR) 
of all events over the T90 and T50 durations by fitting 
their spectra with the spectral functions described in 
section 14.11 These ratios were computed between the 
100-300 keV and the 50-100 keV energy bands from de- 
convolved spectra for direct comparisons with the GRBs 
in the BATSE 4B catalog (jPaciesas et al.l 119990 . Fig- 
ure [5] displays the HR values of the three GBM events 
superposed on the HR diagram of 1973 events with well 
determined durations, HRs and fluences from the 4B cat- 
alog. We notice that for all three events the HRs are 
well within those for the BATSE short GRB population, 
while the durations have a wider spread. For two of the 
events, the Tgo values are in the overlapping area be- 
tween long and short GRBs, while their T50 values are 
well within the short GRB domain. This spread demon- 
strates that the former duration measure is very sensitive 
to 5% background variations, while the latter is more ro- 
bust. Figure [5] also demonstrates that all three events 
are within the short GRB range, albeit, two of them at 
the long tail of the distribution. 




Fig. 5. — Hardness-duration diagrams (Left Panel: Tgo; Right 
Panel: Tso) for 1973 events of the 4B BylTSE catalog (black dots). 
The three GBM short GRBs are overplotted: GRB 090227B (dia- 
mond), GRB 090228 (circle) and GRB 090510 (square). 

We note that a third class of intermediate GRBs has 
been suggested. iHakkila et al.l (|2004[ ). however, ques- 
tioned the interpretation as a physical class, finding, in- 
stead, that the putative class is caused by instrumen- 
tal selection effects. Following the classification scheme 
of Horvath ct al. (2006.), ah three GRBs faU in the short 
class of GRBs. 

4. OBSERVATIONS: II. SPECTRAL ANALYSIS 

We performed simultaneous fits of the TTE data (see 
section [2]) of Nal and BGO detectors selected for each 
of the three short GRBs (see Table [21). We fit the data 
with various spectral functions (see next section) and 
used the Castor C-stat statistic to select the best fit in 
each case. The Castor C-stat differs from the Poisson 
likelihood statistic by an offset which is a constant for a 
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particular dataset; its distribution is similar to a dis- 
tribution for high statistical-content regimes. We chose 
Castor C-stat because of its robustness in the low counts 
regime, which is the case for these 3 GRBs, especially for 
the time-resolved fits of the data. All spectral analyses 
were performed using RMFIT version 3.1 and GBM Re- 
sponse Matrices version 1.8. The detector responses were 
calculated for the best available location of the GRBs, ei- 
ther from Swift or the IPN (see Table [T]). Backgrounds 
were estimated by fitting low-order polynomials to time 
intervals before and after each event. 

4.1. Time- Integrated Spectral Analysis 

Table [S] summarizes the results of our time- integrated 
analysis. For each of the three GRBs we selected intervals 
longer than their T50 durations, as indicated in the Ta- 
ble caption. As in most of the previous spectral analyses 
of short GRBs, we used either the "co mptonized" mode l 
(Compt), or the Band GRB function (jBand et al.lll993[ ). 
The Compt model consists of a power law times an ex- 
ponential cutoff, and has 3 parameters. In our study, 
this cutoff is expressed as -Epcak, namely the lo cation 
of the peak of the vF^ spectrum (jGehrelsl I1997D . We 
also searched for deviations of the spectra from a single 
component by adding a power law (PL) to each of these 
models. 

4.1.1. Identification of an additional component 

For all three GRBs, the combination of Band-f-PL 
provides a significantly better fit than the Band model 
alone. Similarly, for all three GRBs, the combination of 
Compt -l-PL is significantly better than the Compt model 
alone. Table 2] shows the improvement in Castor C- 
stat for the more complex models and the corresponding 
significance of this improvement. All significance val- 
ues reported herein are for comparing nested models, for 
which the probability value can be calculated from the 
difference in Castor C-stat from the distribution for 
the change in the number of model parameters. Test- 
ing Compt-|-PL and Band+PL, one finds that, either the 
two models are tied (GRBs 090228 and 090510) or the 
improvement to the fitting statistic is not large enough to 
justify the more complex Band function (GRB 090227B, 
A Castor C-stat = 3, P = 0.08). 



TABLE 4 

Improvement in Castor C-stat with addition of Power-Law 
component 



Namc3 




M 


odd 




AC-Stat 


S 


gnil 






GRB 090227B 


Co 


Tipt-f PL 


aver C 


□ mpt 


17 


2 


X 


10 


-4 




Ba 


id+PL o\ 


cr Bar 


d 


13 


1 


5 X 


10 


-3 




Ba 


id-l-PL o\ 


cr Cor 


npt-f-PL 


3 





083 






GRB 090228 


Co 


iipt + PL 


jvcr C 


jmpt 


18 


1 


2 X 


10 


-4 




Ba 


id+PL OA 


cr Bar 


d 


18 


1 


2 X 


10 


-4 




Ba 


id-l-PL o\ 


cr Coi 


ipt-l-PL 





1 









GRB 090510 


Co 


iipt + PL 


jvcr C 


jmpt 


25 


3 


7 X 


10 


-6 




Ba 


id-l-PL o\ 


cr Bar 


d 


14 


9 


1 X 


10 


-4 




Ba 


id+PL o\ 


er Cor 


npt + PL 





1 











Figures [S] to [5] show, for each of these GRBs, spec- 
tral fits for Band, Compt-t-PL and Band-|-PL, which 
all provided better fits than Compt (not shown). The 
power-law component appears to overpower the stan- 
dard spectral component below few tens of keV and 



Fig. 6. — Integrated count spectra {left column) and i/Fi, spectra 
{right column) for GRB 090227B with Band row), Compt-|-PL 
(2"<* row) and Band+PL (3'''' row) spectral fits. 

above few MeV, though most of the statistical signif- 
icance is attributable to the fit of the low energy ex- 
cess. The indices of the power laws have similar val- 
ues (~ —1.5) for all events, and they are also simi- 
lar to the value found in the j oint GBM-I-LAT fit t o 
GRB 090510 aheady reported in lAbdo andHI (|2009al ): 
lAckcrmann and al. (2010a) (Table [5]). While there are 
counts in GRBs 090227B and 090510 above the transi- 
tion of the standard function to the additional power law 
at high energies, in GRB 090228, this additional compo- 
nent is only defined by counts below this transition. 

Interestingly, the lag analysis presented in section 13.11 
shows a significant delay of the high energy gamma rays 
above several MeV compared to the low energy gamma 
rays in the case of GRB 090510 and a weak trend in the 
case of GRB 090227B, in agreement with the spectral 
deviation from the Band function identified at high en- 
ergies in these two GRBs. It may be that the positive 
lags identified at high energies in short GRBs are strong 
indicators for the existence of an additional spectral com- 
ponent. Similarly, the weak evidence of time lags mea- 
sured between the 4-20 keV energy band and the other 
ones spreading from 20 keV to 8 MeV for GRBs 090227B 
and 090228 could indicate the spectral deviation from the 
Band function observed at low energy. 

An extra-component abo ve several MeV has be en re- 
ported for the first time in iGonzalez et al.l (|2003[ ) from 
observations of the long GRB 941017 with the Ener- 
getic Gamma Ray Experiment Telescope (EGRET) on- 
board CGRO. Since then, extra components have been 
observed several times in combined PG5'T/GBM-|-LAT 
observations of short and long GRBs, such as 
GRB 090510 (jAbdo and al.l l2009at lAckermann and al.l 
IMOal) and GRB 090902.B (lAbdo and al.ll2009bn . where 
a spectral extension of this high-energy component at 
low energies has also been reported. The possible nature 
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Fig. 7. — Integrated count spectra (left column) and uFi, spectra 
{right column) for GRB 090228 with Band (1''* row), Compt+PL 
(2"'^ row) and Band+PL (3'''' row) spectral fits. 





Fig. 8. — Integrated count spectra [left column) and i/Fi, spectra 
{right column) for GRB 090510 with Band (1=* row), Compt+PL 
(2"'* row) and Band-|-PL (3'"'' row) spectral fits. 

and imphcations of this extra component are discussed 
in section [51 

4.1.2. Companson of the spectral results with the literature 

It has long been known that short GRBs are harder 
than long GRBs, but there has b een a debate about th e 
manner in which they are harder. iPaciesas et al.l ()2003D . 
using a BATSE datatype with coarse temporal (2 s) and 



spectral (16 channels) resolution, found that short GRBs 
were harder in all spectral parameters. When a Compt 
function was used, both the index and i?pcak were harder 
than the ones for the long GRBs; when the Band func- 
tion was used, bot h indic es and ii^ocak were harder. In 
contrast, Ghirland a et al.l (|2004. 2009h) found that the 
higher hardness of short GRBs is entirely due to their 
harder low-energy indices, while their -Epoak distributions 
are similar to long GRBs. 

The -Epoak values of our best time-integrated fits 
(Compt-|-PL) in all events are very high: 2000 ± 100 
keV (GRB 090227B), 720 ± 50 keV (GRB 090228) 
and 3730 ± 270 keV (GRB 090510). The highest two 
values are above any of the other 320 GRBs detected 
with GBM from July 2008 to October 2009, and in 
fact higher than any -Bpoak value reported thus far 
with the possible exception of the Konus/Win d on- 
line catalog of short GRBs |Mazets et a,l.l 1 2004 ■ se e 
also http://www.ioffe.ru/LEA/shortGRBs/Catalog/ 1. 



About 25% of the prompt emission spectra for 86 short 
GRBs with acceptably constrained fit parameters in 
the Konus catalog have -Epeak values, which could be 
above 1 MeV within the error bars, but all are below 
2.5 MeV. Our GBM results show that some bright 
short GRBs have exceptional i^pcak values, beyond 
the range observed in long GRBs. That Epeak val- 
ues for long GRBs end at ~1 MeV is supported by 
the analysis of Sola r Max imum Mission (SMM) data 
by [Harris and S hare' (1998"), who performed a dedicated 
GRB search to address this question. They searched 
for hard-spectrum long GRBs in the 0.8-10 MeV range 
and only 2 were found in the 0.8 - 10 MeV range which 
were not part of the 177 burst sample detected in the 
0.35 - 0.8 MeV range. They concluded that there is 
no population of hard-spectrum long GRBs hidden 
by instrumental effects by SMM, BATSE, and other 
instruments whic h triggered on lo w-ene rgy gamma-rays. 
The analysis of IPaciesas et al.l ()2003D did not reveal 
such extreme values in short GRBs because of the 
limited BATSE/LAD spectral resolution and energy 
coverage (< 2 Me V). An independent analysis of the 
same BATSE data (jGhirlanda et al.l[200l l2009bD found 
a narrow distribution of low i?peak values for all GRBs 
with fluxes below 4 x 10~^ erg.cm~^.s~^, and that short 
GRBs are spectrally similar to the first one or two 
seconds of long GRBs. Our (limited) sample of three 
short GRBs with fluences above 2 x 10^^ erg.cm"^ does 
not support their results. 

Previous analyses of short GRBs have typically 
used a cutoff power law for most fits in prefer- 
ence t o the Band f miction comniqnly used for long 
GRB s (iGhirlanda et al.ll200l I2009bl; IM azets et al]l20Cil 
|http://www .ioffe.ru/LEA/shortGRBs/Catalog/). It has 
not been clear whether the inability to show that a 
Band function is preferred over the cutoff power law in- 
dicates an intrinsic difference between short and long 
GRBs, or merely a lack of s ufficient counts in short 
GRBs (jGhirlanda et al.l l2009al) . This issue is best ad- 
dressed with the broad energy coverage of the GBM de- 
tectors and our selection of a sample that is extremely 
bright, which should enable better determination of the 
high energy power-law index, (3, of the Band function. 
However, we have, as in previous experiments, found it 
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TABLE 5 

Fit parameters of the time-integrated spectra of GRB 090227B, GRB 090228 and GRB 090510 for 4 models: Band function 

(lE. Band), exponential cut-off power law (ie. Compt for Comptonized), Band+Power Law (ie Band+PL) and 
CoMPT-fPowER Law (ie. Compt+PL). The time-interval over which the fits are performed respectively —0.128 s to 0.384 s, 
-0.128 S TO 0.512 s AND 0.500 s to 1.000 S since the trigger time for GRB 090227B, GRB 090228 and GRB 090510. All the 
asymetrical errors are given at the l-cr confident level. The normalization values are computed at 100 keV for Compt 

and at 10 MeV for PL. 



Name 


Model 


Band or Compt 


PL 


Castor C-stat/d.o.f. 






Epeak (keV) 


a or index 




index 




GRB 090227B 


Band 


2116 ± 97 


-0.50 ±0.02 


•^••^^-0.39 


n/a 


699/607 




Compt 


2227 ± 90 


-0.52 ±0.02 


n/a 


n/a 


706/608 




Band-I-PL 


1947l^°« 




_o 44+0.58 
■^•^^-0.80 


-1.51 ±0.05 


686/605 




Compt -l-PL 


1995 ± 96 


-0.36 ±0.05 


n/a 


-1.37 ±0.06 


689/606 


GRB 090228 


Band 


860 ± 50 


-0.59 ±0.03 


o 77+U.b4 
' ' -6.64 


n/a 


813/728 




Compt 


862 ± 52 


-0.59 ±0.03 


n/a 


n/a 


813/729 




Band-I-PL 


723 ± 45 


-0.24 ±0.10 


_4.74l^i4 


-1.64 ±0.03 


795/726 




Compt -l-PL 


722 ± 47 


-0.23 ±0.11 


n/a 


-1 63+"'°^ 


795/727 


GRB 090510 


Band 


4383 ± 290 


-0.75 ±0.02 


-2.80_o 28 


n/a 


911/850 




Compt 


47971^'? 


-0.77 ±0.02 


n/a 


n/a 


922/851 




Band-fPL 


3695 ± 248 


-0.51 ±0.08 


-3.65t«J^ 


-1.38 ±0.04 


897/848 




Compt -l-PL 


3731 ± 265 


-0.58 ±0.07 


n/a 


-1.35 ±0.04 


897/849 


(GBM-I-LAT) 


Band-I-PL 


3936 ± 280 


-0.58 ±0.06 


-2.83 ±0.20 


-1.62 ±0.03 





difBcult to constrain /?. We suggest that this difficulty 
has 3 main causes: (i) Ep^ak is higher in short bursts 
than in long bursts, leaving the measurement of P to be 
made at higher energies where there are fewer counts, 
(ii) P is intrinsically very steep, often resulting in a lack 
of preference for a power-law over an exponential cut-off 
to the spectrum above i?peak, and (iii) the presence of 
an underlying power-law component in addition to the 
standard Band function makes /3 difficult to quantify. 
We conclude that we cannot resolve this question us- 
ing a single spectral component, but instead that an ad- 
ditional power-law component is statistically preferred. 
The Band+PL model is not statistically required, com- 
pared to Compt+PL, but using this model allows us to 
measure /3 or at least to measure upper limits (Table [S]). 
These values confirm the steepness of /?: — 3.44^° ^^ for 
GRB 090227B, < -3.60 for GRB 090228 and < -2.90 
for GRB 090510. 

Finally, the GBM results for the time-integrated low- 
energy ind ex, a, are consistent with prev ious mea- 
surements (jPaciesas et al.l l2003t iGhirlanda et al. 2004^ 

f)9y), and similar to about half of the iMazets et alj 
1041) catalog of short bursts, with a harder than —0.67. 
We conclude here that both a and i^peak contribute to 
the hardness of short GRBs, but the steep high-energy 
spectral-indices indicate a soft spectrum above the break 
energy. 

4.2. Time-Resolved Spectral Analysis 

Dividing the data into fine time intervals reduces the 
number of counts in each spectrum and it is not possi- 
ble to fit models as complex as were used for the time- 
integrated analysis. We selected our time bins to be 
as fine as possible, while still allowing statistically sig- 
nificant spectral analysis by having enough counts dis- 
tributed over the energy channels; the resolution ranges 
from 2 ms to 94 ms. We fit our time-resolved spectra 
only with the Band function since the statistical content 
in small time bins is too low to measure the effect of 
the extra spectral component. Using Band only instead 



of Band+PL artificially makes the high-energy index of 
the Band function softer and i^peak slightly harder. How- 
ever, the effects on the spectral parameters are negligible, 
and in the time bins where Band+PL can be fit, we find 
that the systematic discrepancies between the Band and 
Band+PL fit values are smaller than the statistical er- 
rors. In many cases, we find /3 consistent with — oo and 
provide only la upper-limits. The fits parameters for 
each time interval are presented in Tables [Bl [3 and H] for 
GRBs 090227B, 090228 and 090510, respectively. 

In the past, we have automatically determined time in- 
tervals by accumulating data until a signal-to-noise ratio 
requirement is reached, but this method merges peaks 
and valleys of the light curves, which we strongly wish to 
avoid. In addition, the number of counts required to per- 
form a good fit depends on the spectral shape. Spectra 
spread over a large energy range (such as for hard spec- 
tra) require more counts to have a good fit than spectra 
covering a smaller energy range (soft spectra). We have 
selected statistically significant and variable time bins 
for our time-resolved spectra to account for the various 
pulse structures in the 2 ms light curves. This results in 
time bins containing either peaks or valleys of the light 
curves but not both when possible. To evaluate the relia- 
bility of the spectral fits in these time bins, we performed 
simulations for several low fluence spectra with low and 
high i?poak values as well as with steep /?. The simulated 
sets included backgrounds from the real data; the signal 
was generated using the model function of the real data 
fit. This input function was folded through the GBM 
Nal and BGO responses giving a source count spectrum 
as would be observed with the detectors. Poissonian 
fluctuations were then applied to each energy channel 
of the sum of the source and background for each sim- 
ulated spectrum. Each spectral simulation set included 
ten thousand spectra which were then fit with the in- 
put function model using RMFIT. We found that all fit 
parameters were correctly reconstructed with good con- 
straints and with a systematic bias smaller than the sta- 
tistical errors. 
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TABLE 6 

Parameters of the Band function for the fine- time 

RESOLVED SPECTROSCOPY OF GRB 090227B. AlL THE 

assymetrical errors are given at the l-cr confident level. 
When the high-energy spectral index /3 cannot be 
constrained, the 1-a upper limit is reported. 



TABLE 8 

Parameters of the Band function for the fine- time 
RESOLVED spectroscopy OF GRB 090510. All the 
assymetrical errors are given at the 1-(t confident level. 
When the high-energy spectral index (3 cannot be 
constrained, the 1-a upper limit is reported. 



Tsta 



Tstop 



-0.022 
-0.012 
-0.002 
0.004 
0.010 
0.012 
0.016 
0.018 
0.022 
0.026 
0.028 
0.032 
0.034 
0.038 
0.040 
0.044 
0.048 
0.050 
0.058 
0.076 
0.090 
0.104 
0.112 
0.118 
0.120 
0.126 
0.142 



-0.012 
-0.002 
0.004 
0.010 
0.012 
0.016 
0.018 
0.022 
0.026 
0.028 
0.032 
0.034 
0.038 
0.040 
0.044 
0.048 
0.050 
0.058 
0.076 
0.090 
0.104 
0.112 
0.118 
0.120 
0.126 
0.142 
0.188 



298 + ?i 
— 57 

-169 

^^*"_140 
,+197 



B91 
784 



'--169 
,+208 
-166 
,+292 



09+"-" 

"■"■-'-0.38 

-0.29 ± 0.09 

('■lS_o.l6 
0.02 ± 0.11 



1770if/3 
1387+?« 
2129t«^ 

2482 + «g 

2040t«| 

2916+1°*° 

1711 + 4°^ 
' -307 
1399^291 

^7^0+276 
-197 

741 + "" 
'*^-143 
+ 41 7 

*'°**-126 
1 639+^1° 

«3L\«oi 
189ig 

i°i±r7 



, + 1.82 
-0.81 



1.83^ 
0.11 ± 0.14 
0.71 + 



,+0.67 
-0.32 
0.01 ± 0.16 
0.15 ± 0.15 
0.29 ± 0.29 
0.03 ± 0.18 

°-l'^-0.22 
-0.18 ± 0.12 
-0.71 ± 0.10 
-0.32 ± 0.13 
-0.15 ± 0.19 
-0.51 ± 0.19 
-0.18 ± 0.16 
-0.24 ± 0.12 
-0.47 ± 0.08 
-0.17 ± 0.17 
-0.53 ± 0.08 

02 + °-*= 
°"^-0.31 

-0.52 ± 0.20 
-0.07 ± 0.04 

"■'-''-0.23 

n + 0.58 
-0.37 



; -2.85 
C -4.30 
; -3.16 
C -4.10 
-2.41 
; -3.48 
; -2.88 
; -2.38 
C -4.05 
C -3.91 
; -2.99 
; -2.69 
C -3.45 
C -2.75 



+ 0.30 
0.49 



-3.50 
-2.79 
-2.32 
-2.49 

< -2.92 
-2.37 

< -3.4 
1.96 



+ 0.83 
2.79 
+ 0.71 
1.53 
+ 0.25 
0.35 
+ 0.28 
0.59 



+ 0.29 
0.59 



+ 0.21 
0.31 
2.15 



-0.40^ 



< 1.74 

< -2.36 
-2.00 



+ 0.22 
0.56 



TABLE 7 

Parameters of the Band function for the fine- time 

RESOLVED SPECTROSCOPY OF GRB 090228. AlL THE 

assymetrical errors are given at the 1-a confident level. 
When the high-energy spectral index /3 cannot be 
constrained, the 1-a upper limit is reported. 



Tstart 


Tstop 










el 


-O.OOS 





012 


153 ± 41 


+0 


'^■'-0.43 


-2 


+U.33 
•^■^-0.60 


0.012 





020 


1353 ± 233 


-0 


61 ± 0.08 


< - 


-3.28 


0.020 





022 


l«^±f3'2 


+ 


07 ± 0.22 


< - 


-3.60 


0.022 





024 


1087 ± 237 


+ 


02 ± 0.22 


-2 


46 + 0-3= 
*"'-0.85 


0.024 





026 


„r,,+ 198 
= ■'^-150 


+ 


09 ± 0.26 


-2 


+ 0.47 
°'^-0.88 


0.026 





028 


653 + 1*'= 

"•'■'-lis 


+ 


,,+0.42 
^ -0.33 


-2 


*''~0.62 
.g+0.78 
"''-2.90 


0.028 





032 


738 ± 111 


+0 


08 ± 0.19 


-3 


0.032 





036 


1046 ± 155 


-0 


01 ± 0.16 


< - 


-3.97 


0.036 





042 


711 ± 168 


-0 


+ 0.27 
^'■^-0.18 


-2 


94 ± 2.68 


0.042 





048 


623 + 1"^ 
t>^3_89 


-0 


21 ± 0.17 


< - 


-2.65 


0.048 





054 


575 ± 88 


+ 


05 ± 0.24 


< - 


-3.08 


0.054 





060 


479+^= 


-0 


03 ± 0.22 


< - 


-3.45 


0.060 





066 


597 ± 76 


+ 


31 ± 0.23 


< - 


-3.93 


0.066 





072 


728li;= 


-0 


03 ± 0.18 


< - 


-2.99 


0.072 





078 


700 ± 119 


+0 


02 ± 0.20 


< - 


-2.56 


0.078 





082 


451 + lf 


+ 


63 ± 0.48 


-3 


n,+0.61 
"'^-2.73 


0.082 





092 


650 ± 94 


+ 


09 ± 0.20 


< - 


-3.01 


0.092 





128 


307 ± 35 


-0 


14 ± 0.21 


< - 


-3.19 


0.128 





170 


42 ± 14 


+ 


g2 + 2.41 
■^^-1.67 


-2 


92+0.45 
■'^-0.83 


0.170 





212 


23 ± 24 


-1 


= , + 0.00 
^^-0.46 


< - 


-2.68 



4.2.1. -E'pcak distribution and evolution 

The distribution of the time-resolved £^pcak for the 
three GRBs ranges over a very wide range of values (from 



Tstop 



0.5 

0.52 

0.53 

0.534 

0.538 

0.540 

0.542 

0.546 

0.548 

0.55 

0.554 

0.56 

0.566 

0.568 

0.570 

0.576 

0.588 

0.592 

0.596 

0.600 

0.604 

0.608 

0.612 

0.616 

0.62 

0.632 

0.642 

0.652 

0.658 

0.664 

0.67 

0.678 

0.684 

0.692 

0.710 

0.714 

0.726 

0.736 

0.756 



0.52 

0.53 

0.534 

0.538 

0.540 

0.542 

0.546 

0.548 

0.55 

0.554 

0.56 

0.566 

0.568 

0.570 

0.576 

0.588 

0.592 

0.596 

0.600 

0.604 

0.608 

0.612 

0.616 

0.62 

0.632 

0.642 

0.652 

0.658 

0.664 

0.67 

0.678 

0.684 

0.692 

0.710 

0.714 

0.726 

0.736 

0.756 

0.850 



SSIJ 
418"^ 



,0.0+1190 
^■^^■^-607 

+ 2030 
■=i^i3'^i_1200 

2460 + 1"^" 
^^'^"-669 
,+294 
-202 
+ 160 
'-112 
^.-.+1493 

2662"'"^^^ 
1382l™° 
1010±f,3 
2397+^36 

2233 + 3|| 

,.„.+598 
3*°"-484 
394n+''25 
3^*"-540 
1146 +==3 
^^*'0-344 

5181 + 112° 
■''-''^-4210 

3392 ± 855 

2606 + 1°*'° 
20UO_ J290 
1013 + 3560 
1U13_493 

+ 1630 
^iUl_g27 
,„S7+1420 

7641 + 3060 

2443 + 1430 

+ 6740 
*°°^-1730 
4699 + 1""° 

+ 1320 
'>^l'J-989 
3369 ± 485 
,„, . + 1450 
3 = 14-889 

3576 + 1*13 
3="'-1064 
- + 1520 
■1170 
, + 593 
"-491 
, + 1526 
-1130 
2771 ± 836 

355 + 216 
"■"^-131 
1 + 219 

-53 



«+U.77 

-0.43 
+ 1.65 
-0.63 
fO.35 
-0.25 



-0.68^ 
-0.18 + 
-0.74j 
-0.61 ± 0.17 
-0.80 ± 0.14 
-0.70 ± 0.17 
-0.36 ± 0.24 

-0 47+0-93 
O-«'-0.56 

-0 57+0-4= 
O-"'-0.28 

-0.73 ± 0.13 
-0.48 ± 0.13 
-0.11 ± 0.27 

+ ,2 + 0-43 
+ 0-32_o.62 
^ ^^ + 0.32 
-0.66_o 23 
-0.62 ± 0.16 

-0.27 ± 0.15 

o.5oj; 

-0.19 ± 0.28 

+0.06l0-6t 
-0.91 ± 0.25 

-0 71+1-26 
o- '1-0.60 

-0 16+0-6= 
•J -^0-0.32 
-0 14+0-''2 
-0.55 
-0.72 ± 0.32 



rn+0- = l 
'"-0.37 



< -2.26 

< -2.18 

< -2.56 

< -3.18 
-2.51±°-=^ 

< -2.94 

< -3.09 

< -2.80 

< -2.01 

< -2.31 

< -3.12 

< -3.57 

< -2.81 

< -2.54 

< -3.10 

< -4.23 

< -3.46 

< -3.69 

< -3.17 

-1 73 + °-44 
^- '^-0.83 



,+0.36 
-0.21 



4445^ 
19S2j 
6646^ 



159^ 



-0.63_ 

-0.70 ± 0.19 

+ 1.34 ± 1.20 

+ 14+1-20 
+ "--'*-0.34 

+ 46+l-°l 
+ U-'>6_o 3g 

-0.39 ± 0.16 
+ 0.05 ± 0.27 
-0.44 ± 0.22 
-0.30 ± 0.32 
-0.32 ± 0.27 
+ 0.97^ 
-0.21 ± 0.25 
-0.41 ± 0.20 
-0 37+0.43 
0-3'-0.29 
-0 23+0-64 



-2.66 
-1.71 

< -2.59 

< -1.77 

< -2.28 



+ 0.65 
1.32 
+ 0.23 
0.90 



-2.54 



+ 0.37 
0.84 



+ 0.45 
0.84 



7+1.73 
-0.76 



-2.3 
-2.60 

< -3.61 

< -3.65 

-2 92 + °-61 
''■"2_5 11 

-2 61 + ° =* 
'^■"'--1.19 

-2 69+°-=0 
^■'^^-2.96 

< -2.87 

< -2.91 
-2.64±°- = ; 
-1.74 ± 0.17 
-1.51 ± 0.07 



100 ± 30 keV to 3.2tJ'? MeV for GRB 090227B, from 
23 ± 24 keV to 1.5 ± 0.3 MeV for GRB 090228 and from 



I59I53'' keV to 7.6Z2 MeV for GRB 090510), indicatmg 
a strong spectral evolution in all events (see top panel of 
figure[5]). Such a strong spectral evolution is very unusual 
compared to what one observes in long bursts, where the 
-E'pcak distribution does not extend to the high energies 
found here (Ford et al. 1995). 

Panels c and d of figures [U [2] and [3] show the evolution 
of ii'pcak compared to the count-rate light-curves in two 
energy bands, 8 - 200 keV and 1-38 MeV, with the 
same time bins used to compute ii'pcak. For all three 
GRBs, the i?poak values follow a global soft-hard-soft 
evolution as also seen in the high temporal resolution 
light-curves presented in section [31 We notice that E'pcak 
mostly tracks the light curves, approximately correlating 
with the shape of the count-rate variations but not al- 
ways with the intensity. We quantified the E'pcak ~ inten- 
sity correlation by computing the Spearman Rank Order 
Correlation Coeffient in the two energy ranges described 
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Fig. 9. — Top: Distribution of the peak energy, iiJVioaki of the 
Band spectral function for the same 3 GRBs. Bottom: Distribution 
of the l-CT lower limits of the low energy index a of the Band 
spectral function for GRB 090227B, GRB 090228 and ORB 090510. 
The vertical lines indicate the electron synchrotron emission limits 
for fast cooling (index of -3/2) and slow cooling (index of -2/3). 

above (See Table We find that the correlation is most 
significant in the highest energy band for all three GRBs. 
Even with this small sample, the independent correlation 
measurements increase the robustness of our results. 



TABLE 9 

Spearman Rank Order Correlation Coefficients (p) and 

ASSOCIATED PROBABILITY COMPUTED BETWEEN THE RAW COUNTS 
LIGHT CURVES IN THE 8-200 KEV AND 1-38 MeV ENERGY BANDS, 
AND THE i?poak VALUES PLOTTED AT THE BOTTOM OF 
FIGURES [1] TO [3] 



Energy 


Dands : 


8 - 


200 kcV 


1 - 


38 McV 


Name 


Nu 


iibcr of data points 


P 


probability 


P 


probability 


GRB 090227B 




27 


0.354 


0.076 


0.452 


0.021 


GRB 090228 




20 


D.49S 


0.030 


0.666 


0.002 


GRB 090510 




39 


-0.014 


0.935 


0.371 


0.022 



The ji^Dcak-intensity correlation is a well known charac- 
teristic of long GRBs ()Ford et al.ll 19951 iGolenetskii et al.l 
[1983), wher e frequently the p eak in i?peak leads the peak 
in intensity (jFord et al.lll995| ). This pattern seems weak 
or absent in our short events, possibly because the low 
statistics in such short time intervals make it necessary 
to combine too many bins to see this effect. To study 
the simultaneous changes in Speak values and intensities 
of these events we computed and displayed the deriva- 
tives of their values (we cannot use cross-correlation tech- 
niques as we are limited by the small number of data 
points). Figure [TU] shows a scatter plot of these deriva- 
tives for all three bursts. Points in the +/+ and -/- 
quadrants indicate a correlated increase or decrease of 
the two data sets, respectively. A positive correlation is 



1200 — 
1000 — 



-20000 -10000 10000 20000 30000 40000 50000 

A Flux (counts/s) 

Fig. 10. — Top: Variation of the peak energy -Epoak versus 
the count rate variation for two energy bands 8-200 keV (stars) 
and 1-38 MeV (triangles) for GRB 090227B, GRB 090228 and 
GRB 090510 (red, blue, and green, respectively). Bottom: zoom 
in on GRB 090228 only (8-200 keV; blue stars; 1-38 MeV: red 
triangles) . 

evident, especially in the 1-38 MeV light curves, indi- 
cating that increasing (decreasing) flux correlates with 
increasing (decreasing) iJpcak- 

4.2.2. Low and high-energy power-law indices 

The bottom panel of Figure IH] shows a histogram of 
the Icr lower-limits of the low-energy index a. The a 
values are exceptionally hard, increasing the discrepancy 
with the predictions of synchrotron emission even if they 
don't rule out this emission mechanism as the main pro- 
cess to explain the prompt emission, as we will discuss 
in section [S] The time-resolved values of a are even 
harder than the time-integrated ones (also described in 
section l4.1.2|) and most of the values exceed — |, the elec- 
tron slow-cooling limit. 

This very hard distribution f or both a and ij^pea k 
confirms the results presented in iPaciesas et al.l ()2003() , 
which describes the short GRBs as harder with both 
higher i?pcak values and harder a, and does not sup- 
port the explanation offered by IGhirland a et al. (2003, 
l2009b). 

It is also remarkable that the time-resolved (3 val- 
ues (Tables [H] to [5]) are mostly very soft according 
to their Icr upper limits: from < —4.30 to < —1.74 
for GRB 090227B, from < -3.97 to < -2.11 for 
GRB 090228 and from < -4.23 to < -1.29 for 
GRB 090510. The variations of this index between 
the hardest well-constrained values and the softest la 
upper limits also indicate strong evolution: from < 



-4.30 to -1.96 



-0.21 
-1-0.31 



for GRB 090227B, from < -3.97 
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to -2.39i^;^g for GRB 090228 and from < -4.23 to 
-1.73+°:^^ for GRB 090510. 

5. INTERPRETATION 

The prompt emission from short and fong GRBs shows 
remarkably similar properties. The main difference ap- 
pears in the hardness- duration diagram : short GRBs 
are harder (jKou vel iotou et al.l [l993). It is therefore 
usually believed that short and long gamma-ray bursts 
are produced by different progenitors leading to the 
same succession of events: formation of a compact 
central engine and relativistic ejection in the form of 
collimated jets (e.g. Nakar 2007, for a review). In this 
picture, the prompt emission is produced by the same 
mechanism for all GRBs and the observed differences 
between short and long bursts have to be explained 
by different initial conditions (energetics and lifetime 
of the central engine, source environment, etc.). Due 
to the short variability timescale observed in GRBs, 
the prompt em i ssion has very likely an internal origin 
(|Sari and PiranI 119971 ). which means that it is not 
related to the deceleration of the jet by the ambient 
medium. This scenario favors a greater similarity in 
the prompt properties of short and long GRBs than for 
their afterglows. There are then three possible energy 
reservoirs from which the prompt radiated energy can be 
extracted, in a region located above the photosphere of 
the relativistic outflow and below its deceleration radius: 
(i) the internal (thermal) energy, that can be extracted 
at the photosphere (e.g. Pe'er 2008: B eloborodov|[200l 
for recent discussions of this mechanism) ; (ii) the kinetic 
energy, that can be extrac ted via (internal) shock waves 
(|Rees and Meszaroill994D : (iii) the magnetic energy - 
if the outflow is strongly magnetized - which can be 
extracted via magnetic reconnectio n (e.g. iSpruit et al.l 
I200H iLvutikov and BlandfordI 120031 ). In the latter case, 
the physics of the magnetic reconnection in relativistic 
outflows is still far from being understood and it is 
difficult to predict precisely the spectrum of the expected 
radiation. For this reason, we will not discuss it further. 
In addition to the mode of energy extraction, one also 
needs to understand the distribution of the emitting 
particles, which constrains acceleration mechanisms, 
and the nature of the dominant radiative processes. 



5.1. Why are short GRBs harder? 

The results presented here offer new clues to identify 
which mechanism is at work during the prompt emission 
of short GRBs. In the three bright short GRBs studied 
here, the broad spectral range of the GBM instrument 
on board FGST allows us to demonstrate clearly that 
the origin of their hardness is mainly due to very high 
peak energies, well above the usual range observed in 
long GRBs. This is naturally ex pected in internal shocks 
(|Daigne and Mochkovitch|[l998() . Precisely, the observed 
isotropic peak energy, E ^^^^^, associated with a collision 
scales as ()Barraud et al.ll2005l ) 

,1/2 2 pi/2 

i^^bs'x^0(^) (1) 



where E and T are the isotropic equivalent kinetic energy 
flux and mean Lorentz factor in the outflow and r is the 
variability timescale during the relativistic ejection. The 
contrast k is deflned as the ratio between the maximum 
and the minimum Lorentz factor in the outflow a nd the 
function 4>{k) is steadily increasing (see Barra ud et al.l 
|2005[) . The microphysical parameters Ce and cb are the 
fraction of the energy dissipated in internal shocks that 
are injected respectively in relativistic electrons and in 
the magnetic field, and ( is the fraction of electrons 
that are accelerated. The isotropic equivalent gamma- 
ray luminosities in the small sample of short bursts that 
hav e a measured redshift are simi lar as for long GRB s 
(see lNakaii pOOTi ) and Fig. 3(b) in lZhang et all ((20091) ). 
Therefore, one can assume similar kinetic energy fluxes 
E in short and long GRBs. On the other hand, we show 
here that the lightcurves of the three short GRBs imply 
that the variability timescales are contracted, compared 
to long GRBs, during the relativistic ejection, which is 
most probably caused by the differences in the central en- 
gine. If all other properties are flxed (distribution of the 
Lorentz factor, microphysical parameters), this leads to 
higher peak energies as E^^^^^ oc and gives a natural 
explanation for the observed hardness-duration relation 
in GRBs. This relation should only be a general ten- 
dency, since variations in the Lorentz factor (F, k) or 
the microphysical parameters (^, Ce, ee) from one burst 
to another will lead to some dispersion according to the 
equation above. 

More generally, for many possible radiative processes, 
the peak energy in the comoving frame should scale as 
-^p^com P*^* where and e* are the comoving den- 
sity and specific energy density in the emitting shocked 
region. This leads to the following general expression of 
the observed peak energy (Barraud et al. 2005) : 

4%s°^*-.(^)f55^' (2) 

where ^xy is now a function depending on the two ex- 
ponents X and y. Synchrotron radiation with constant 
microphysics parameters corresponds to a; = 1/2 and 
y — 5/2, which leads to equation [T] above. Inverse 
Compton scattering in Thomson regime corresponds to 
X = 1/2 and y = 9/2, which leads to the same depen- 
dance of the peak energy on the duration. The fact that 
the observed peak energy increases in internal shocks 
when the timescale decreases remains valid as long as 
a; > 0, e.g. as long as the peak energy in the comoving 
frame decreases with the density, which seems a reason- 
able assumption. This is for instance the case for the 
jitter radiation (Medve dev et al.ll2009D if the correlation 
scale of the magnetic field increases when the density de- 
creases. While the hardness-duration relation seems a 
robust feature of the internal shock model, the slope of 
the predicted relation depends, however, on the details of 
the dominant radiative process and/or the acceleration 
mechanism in mildly relative internal shocks, as well as 
the dispersion due to the distribution of the other pa- 
rameters (Lorentz factor, kinetic energy flux, etc). 

In the internal shock model, it is also expected 
that the peak energy should track the lightcurve. If 
a given pulse is associated with the propagation of 
a shock wave within the outflow, the magnetic fleld 
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decreases when the radius increases, resuhing in a decay 
of the peak energy in the decay p hase of the pulse 
paigne and MochkovitchI [19981 [200l . 

In photospheric models, the observed peak energy is 
related to the temperature of the fireball at the photo- 
spheric radius and should therefore scale as 

ijTff ^ . ]° , (3) 

p,obs £'5/12 ' ^ ' 

where rg is the radius where the initial relativistic 
ejection takes place. It appears then less natural to 
expect short GRBs to have higher peak energies in this 
scenario, as the shorter timescales in short GRBs would 
imply smaller radii tq. Higher peak energies would be 
obtained if - due to differences between the central 
engines - outflows in short GRBs have systematically 
higher Lorentz factors. As only lower lim its on F have 
been obtained so far for FGST GRBs (jAbdo and al.l 
[2009b. a, c), this is still a possibility. However, as these 
lower limits are already quite large, this could lead to a 
severe constraint for models of jet acceleration in central 
engines of short GRBs. 

5.2. On the low energy power law index, a 

The simple scenario where the prompt emission in 
(short and long) GRBs is due to emission of relativistic 
electrons accelerated in internal shocks within a rela- 
tivistic outflow, offers a consistent framework to explain 
many observed features: high variability in lightcurves, 
spectral evolution, time lags, etc. In addition, compared 
to the photospheric model, it seems to give a natural 
explanation for the observed hardness-duration relation, 
which is a fundamental observation in the comparison 
between the two classes of bursts. In the internal shock 
model, short GRBs have higher peak energies because 
of shorter variability timescales, which seems consistent 
with the GBM observations discussed in this paper. 
Most of these properties are due to the dynamics of 
internal shocks and are therefore robust predictions, 
found for different assumptions regarding the dominant 
radiative process and the microphysics in mildly rel- 
ativistic shocks. The two most discussed possibilities 
are synchrotron radiation or synchrotron self-Compton 
(SSC). The fact that strong components in the GeV 
range, dominating the total energy release, are not a 
common feature in FGS T bursts favors the synchrotron 
scenario for long bursts (jBosniak et al.|[2009l : iPiran et"al] 
I2009D . There is, however, a remaining worrying problem 
that appears once again in the current analysis: the 
observed low-energy photon index, a, is too steep 
to be easily explained by synchrotron emission. The 
expected slope in the fast-cooling regime is a = —3/2 
(jGhisellini et al.l |2000[) . Inverse Compton scattering in 
the Klein-Nishina regime can aff'ect the electron cooling 
and steepen the sync hrotron sp ectrum up to a = — 1 
(IDerishev et al.l 120011 : [Bosnjak et al.l 120091: INakar et all 
120091) . Svnchrotron radiation in the slow-cooling regime 
allows to reach a = —2/3 but leads to a problem with 
the radiative efficiency. In the present sample of three 
bright short bursts, the measured values of a are almost 
always above this limit a = —2/3, which is another 



differenc e with long GRBs s howing a mean value close to 
a = —1 (jPreece et al.ll200dl) . These observed values of a 
are really challenging for the synchrotro n process. Jitter 
radiation can produce such flat slopes ([Medvedev et al.l 
[200i). It remains to be confirmed, however, that it is a 
viable proc ess for physical condit ions expected in inter- 
nal shocks (jMedvedev et al.ll2009[ ). The SSC mechanism 
can also lead to steeper values of a but it should then be 
understood why the additional component in the LAT 
range are not brighter. Photospheric models can also 
show very steep low-energy spectral slopes but on the 
other hand, they do not propose natural explanations 
for the hardness-duration relation and more generally 
the observed spectral evolution. 



5.3. On the high energy power law index, /3 

The broad spectral range of GBM in principle allows 
investigation of the values of the high-energy slope (3 
but the observed high values of the peak energies makes 
this measurement difficult in the three short GRBs 
presented here. In many cases, a power-law at high 
energies is not statistically preferred to an exponential 
cutoff. The time-resolved spectral analysis presented 
here allows us, however, to constrain the value of /3 in 
most time bins. The measured values of /? are almost 
always compatible with /? < — 2 which shows that the 
main component is recovered, without any need for an 
additional break at higher energies. If the dominant 
process is synchrotron radiation from shock-accelerated 
electrons, such measurements can be translated into 
constraints on acceleration theory in mildlj[3 relativistic 
shocks. This regime is not well covered by present 
simulations which usually focus on the ultra-re l ativisti c 
regimes expected in afterglows (jSpitkovskvl l2008al ). 
The fact that f3 is variable within a GRB would be 
an indication that there is no universal slope for the 
distribution of shock-accelerated electrons in the mildly 
relativistic regime. In addition, the observation that /3 
is very steep would also indicate that the distribution of 
shock-accelerated electrons is very peaked at low Lorentz 
factors, as the predicted slope —p ~ 2 {(3 + 1) lies in the 
typical range —6 to —2 for /3 = —4 to —2. This could 
also simply indicate that the distribution of electrons 
is more complex than a s ingle power-law ([Spitkovskvl 
I2008bt [Martins et aI[|2009D . The true distribution of 
electrons could however be masked by other effects that 
could steepen [3, especially if the real spectral shape at 
high energies for the main component is a power-law 
of slope (3 above -Epoak followed by a cutoff above 
Ecut > -Epeak- The mcasurcd slope could appear as 
steeper than f3 if i^cut is not at too high energy. Several 
physical processes could be responsible for such a cutoff: 
for instance the existence of a maximum Lorentz factor 
for shock accelerated electrons or 77 annihilation above 
the threshold for pair production. The resulting spectral 
shape could be rather complex. It is beyond the scope 
of this paper to model such processes in detail but the 
broad spectral range of GBM, especially when combined 

Unlike the external shock, which is a relativistic collision be- 
tween the relativistic outflow and the resting ambient medium, in- 
ternal shocks occur between density fronts with mildly-relativistic 
relative-velocities within the outflow. 
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with the LAT, clearly offers a new opportunity to make 
progress in this field in the future. 



5.4. On the additional spectral component 

In addition to the main component which is well fit- 
ted by the Band GRB function, there is evidence of 
an additional power-law component, dominant at low 
energy, and - at least in the case of GRB 090510 
and GRB 090227B - an addition al component at high 
energy which could be related (jAbdo and al.l l2009al : 
lAckermann and aI]|2010aD . This suggests that the ra- 
diative processes at work are more complicated than was 
predicted in the simplest model, possibly mixing several 
distinct components with varying relative weights. 

The difficulty in identifying the mechanism responsible 
for the main spectral component that peaks in the GBM 
energy range makes the physical modeling of this addi- 
tional component even more challenging. In scenarios 
where the soft gamma-rays are produced at the photo- 
sphere, the additional component could be the signature 
of synchrotron radiation from fas t-cooling electron s ac- 
celerated in internal shocks (e.g. iRyde et 311120101 ). It 
implies, however, that the peak /ifm of the synchrotron 
radiation is at high energy (LAT range), which leads to 
extreme constraints on the magnetic field B and the elec- 
tron Lorentz factor Fm : 



BVi 



6 X 10' 



—) 

1000 / 



^^m,obs 

IGeV 



G (4) 



On the other hand, if internal shocks are responsible for 
the soft gamma-ray emission, several mechanisms can 
produce additional components at low and/or at high 
energy. 

In addition to the leptonic synchrotron component, a 
weak in verse Compton component can be present at high 
energy (iBosniak et al.|[20 09). as well as a hadronic com- 
ponent ( Asano et al.ll2009f ). The lack of identification 
of a clear cutoff in the high-energy spectrum of GRBs 
detected by the LAT makes it difficult to evaluate if a 
strong 77 annihilation occurs at high en ergies even if 
recent observations indicate such a cutoff ()Abdo and al.l 
I2009bl: lAckermann and al.l [2QTQbl ). If this is the case, 
the new component at low energies could be associ- 
ated with the radiation of the produced electron-positron 
pairs and would then be related to the high-energy 
component. Alternatively, even in the internal shock 
model, photospheric emission is usually also expected, 
and could be superimposed on the synchrotron radia- 
tion (jPaigne and MochkovitchI 120021) , which can again 
lead to an additional component in the spectrum, prefer- 
entially in the GBM range or below. Note that in the lat- 
ter case, different central engines and acceleration mech- 
anisms for the outflow could lead to different intensities 
of the photospheric component in short and hard GRBs. 
We cannot assess here the viability of these models, but 
our analysis does illustrate that introducing more real- 
istic and complex physics is probably now necessary to 
model the spectral shape observed in GRBs. Finally, 
one should keep in mind that our study is based on three 
short bursts only, selected by their high brightness and 
that they are not necessarily fully representative of the 
short GRB population. 



6. CONCLUSIONS 

With GBM, it has been possible for the first time to 
perform time resolved spectroscopy of short GRBs with 
a resolution down to 2 ms. The time-integrated spectra 
of the three GRBs analysed here are best fit with the 
sum of a power law with an exponential cutoff and an 
additional power law with indices similar for the three 
GRBs, clustering around —1.5. This additional power 
law indicates a deviation from the standard Band func- 
tion usually used to fit GRB prompt emission spectra in 
the keV-MeV energy range and can be interpreted as an 
additional component. We find that this extra compo- 
nent overpowers the Band function at low energies (be- 
low a few tens of keV), and at high energies (above a 
few MeV), and could be a combination of various mech- 
anisms as explained in Section [SJ How the low and high 
energy excesses are related together is still not clear yet. 

The spectral parameters obtained from the fits of the 
integrated spectra are much harder than usually observed 
in long GRBs for both the low energy index a (> —0.43) 
and the ii^poak (up to about 4 MeV) values. When a 
Band-f-PL model is used, the high-energy spectral index 
below —2.90 indicates a soft spectrum above the break 
energy. The time-resolved and time-integrated low en- 
ergy indices nearly always violate both the electron slow 
and fast cooling predicted by the synchrotron models and 
additional emission processes may be required to explain 
the data. 

Based in the observed pulse structures, we posit that 
the light curves of short GRBs are similar to the ones 
of long GRBs but contracted in time. Future observa- 
tions that probe even shorter timescales in both types of 
GRBs will test this hypothesis. The fine-time resolved 
spectroscopy shows that the short GRB -Epcak values are 
stretched up to higher energies compared to long burst 
ones. The time-resolved i^poak values can vary from few 
tens of kcV up to several MeV in a time scale as short 
as a few hundreds of ms. i?pcak mostly tracks the light 
curve evolution in a way that is similar to long bursts. A 
convincing hardness-intensity correlation has been mea- 
sured in these three GRBs between -Bpeak and the light 
curves above 1 MeV. A more general correlation between 
the signs of the derivatives of the -Epcak curves and the 
light curves seems to be present for these three GRBs. 
These results favor the scenario where the same physical 
mechanism is at work in the prompt phase of short and 
long GRBs. The fact that short GRBs are harder would 
then be mainly due to their shorter timescales, proba- 
bly a signature of different central engines. The internal 
shock model naturally leads to the expected hardness- 
duration relation for a broad range of radiative processes 
whereas such a relation appears less natural for other 
energy extraction mechanisms. 
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